Charmed and Bottom Baryons from Lattice NRQCD 



Nilmani Mathur"^'', Randy Lewis^ and R. M. Woloshyn'^ 
"■TRIUMF, 4004 Wesbrook Mall, Vancouver, British Columbia, Canada V6T 2A3 
^Department of Physics and Astronomy, University of Kentucky, Lexington, KY 40506-0055 
'^Department of Physics, University of Regina, Regina, Saskatchewan, Canada S4S 0A2 



The mass spectrum of charmed and bottom baryons has been computed on anisotropic lattices 
using quenched lattice nonrelativistic QCD. Masses are extracted by using mass splittings which are 
more accurate than masses obtained directly by using the nonrelativistic mass-energy relation. Of 
particular interest are the mass splittings between spin-1/2 and spin-3/2 heavy baryons, and we find 
that these color hyperfine effects are not suppressed in the baryon sector although they are known to 
be suppressed in the meson sector. Results are compared with those obtained in a previous NRQCD 
calculation and with those obtained from a Dirac- Wilson action of the D234 type. 
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I. INTRODUCTION 

A comprehensive knowledge about the mass spectrum 
and spin sphttings of heavy baryons is important for our 
understanding of quantum chromodynamics. However, 
except for singly heavy charmed baryons and only one 
singly heavy bottom baryon (Ah), most of the heavy 
baryon masses have not yet been measured experimen- 
tally On the theoretical side there are many results 
on heavy baryon masses from different models including, 
for example, a number of quark model variations^, ^, Q. 

Using lattice QCD, substantial work has been done in 
the heavy meson sector. However, so far only very few 
results have been reported for heavy baryons[^ ^l; 0; Hi 
and there is only one work Q where heavy (bottom) 
quarks are treated nonrelativistically. A further study of 
charmed and bottom heavy baryons on the lattice using 
nonrelativistic QCD (NRQCD) therefore seems worth- 
while. 

Extraction of the experimentally observed mass split- 
tings between vector and pseudoscalar mesons remains 
a challenging problem in lattice QCD; quenched calcula- 
tions have so far been unable to extract the observed mass 
splittings [pi , and unquenched studies have not resolved 
the issue Therefore, it is natural to ask whether lat- 
tice results for baryon mass splittings also exhibit similar 
suppression compared to experiment. 

Empirically, spin splittings in baryons are smaller than 
those in the meson sector. Moreover, in a lattice simu- 
lation the correlators for baryons, particularly for spin 
3/2 states, are noisier than those for mesons, and thus, 
by using lattice QCD it is comparatively difficult to ex- 
tract a reliable mass spectrum for heavy baryons. In this 
work we report on the charmed and bottom baryon mass 
spectrum and mass splittings by using a nonrelativistic 
heavy quark action and an improved light quark action 
on anisotropic lattices. 

In section II we summarize different charmed and bot- 
tom heavy baryons with their relevant quantum numbers 
and discuss our choice for interpolating fields. Section III 



presents numerical simulation details. For heavy quarks, 
we use the nonrelativistic action from Rcf. while a 
tadpole improved gauge action and an improved Dirac- 
Wilson action of the D234 type are used for light 
quarks. Since these actions were previously detailed else- 
where H, we will describe them only in an appendix. 
The calculations are done on two different anisotropic 
lattices with the same gauge configurations as were used 
in Ref. § at /3 = 2.1 and /3 = 2.3. 

In section IV we present our results. Masses are cal- 
culated using two methods; the first uses the standard 
NRQCD relation between mass and energy while the sec- 
ond employs mass splittings to calculate masses. As mass 
splittings can be estimated more accurately than masses, 
errors in the second method are smaller than those ob- 
tained from the first one. The overall systematic uncer- 
tainty is estimated by including scale uncertainty, un- 
certainty due to the choice of a time window for fitting 
correlation functions, error due to extrapolation to the 
physical light quark masses, uncertainty in fixing charm 
and bottom masses and uncertainty from our determina- 
tion of the lattice anisotropy. 

Spin splittings are discussed in section V. From our 
results, along with other published results, we conclude 
that the suppression of mass splittings is not present in 
the baryon sector in the same way as it is in the me- 
son sector. Over the whole mass range where data are 
available, quenched lattice QCD simulations yield mass 
differences between spin 3/2 and spin 1/2 baryons which 
are comparable to or larger than experimental values. 

II. CHARMED AND BOTTOM BARYONS 

Singly and doubly charmed and bottom baryons are 
summarized in Tables I and II respectively. Table II 
also includes doubly heavy states containing two differ- 
ent heavy quarks (charmed and bottom quarks together). 
Quark content, as well as the spin-parity J^, the isospin 
/, and si which identifies the total spin of the light quarks 
(also spin-flavor symmetry: s; = is symmetric while 
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TABLE I: Summary of singly heavy baryons, showing valence 
quark content {q = u,d and Q = c,b), spin-parity, isospin and 
mass (in GeV). The quantity si is the total spin of the light 
quark pair. The experimental values are from Ref. 
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TABLE II: Summary of doubly heavy baryons, showin; 
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s/ = 1 is antisymmetric) are shown. Notice that masses 
for many singly heavy states are not measured yet and 
there are no data at all on masses for doubly heavy states. 

To project out heavy baryon states we use the same 
interpolating operators as were used in Ref. |^. For S- 
like baryons we choose 

E : e'^^^feJC'TsQfclgc (1) 

where g is a light quark field and Q is a heavy quark 
field. Here a, 5, c are color indices whereas Dirac indices 
have been suppressed. For Sg, q is u or d and for Oq, q 
is s. For doubly heavy E-like baryons with equal heavy 
masses, we interchange the role of light and heavy fields, 
i.e., to get Sqq, we change q ^ Q and Q — > u or rf. 
Similarly, for ilgg, the change is g — > Q, Q — > s. 

The Sq is E-like but it contains two different light 
flavors so it is considered separately as 

S' : ±{e^''-[q'TCj,Q,]q, + e'^'^[q'^CjM,} A'^) 

with q — u or d and q' = s. 

The A-like baryons involve three distinct flavors. A 
simple choice is the heavy lambda : 

A : e^'^iq^C-fM,. (3) 

where for Aq, q — u, q' — d, and for Sq, q = u, q' ~ s. 
A more symmetrical choice would be the octet lambda 

-[qaC^5Qb]qcy (4) 

with the same flavor assignment as for the heavy lambda. 
One can use either of these A states as they give con- 
sistent results^. We choose the octet-lambda (Aq) for 
this work. For spin 3/2 baryons we choose the following 
interpolating field: 

E* : e'^'^iq^CjMc (5) 



where for Eg, q = q' is u or d and for fig, q — q' is 
s. To get 5q, one needs to consider q = u or d and 
q' = s. Similarly, to get the doubly heavy states with 
equal heavy masses one needs to interchange the role of 
light and heavy fields. For example, to get Eqq, one 
needs to change q,q' ^ Q and Q — s- u or c?, whereas Q,*q 
requires q,q' — > Q and Q ^ s. 

The operator in Eq. (||) has both spin 1/2 and spin 3/2 
states. At zero momentum the corresponding correlation 
function can be written as [ p3[ : 

Cij{t) ^ ((5y - ^7j7j)C3/2(i) + ^liljCi/2{t), (6) 

where i,j's are spatial Lorentz indices and C3/2(i/2) 
are the spin projections for spin 3/2 (1/2) states. By 
choosing different Lorentz components the spin 3/2 part, 
^3/2(^)7 is extracted and used to calculate the mass of 
the spin 3/2 baryons. 

Operators for baryons with two unlike heavy flavors 
may be constructed from the above interpolating opera- 
tors by interchanging the role of heavy and light fields. 
For example, Sgg, and ^qqi can be obtained from 
Eq. (5) by letting g, q' — > Q, Q' and Q q with q — u or 
d and q = s respectively. For Sqq, and ^'qqi we use the 
symmetrical form again, as given by Eq. (3), making the 
same replacements i.e., q, q' — + Q,Q' and Q ~^ q with 
q = u or d and q = s respectively. Finally, Sqq' and 
^QQ' ^re the doubly heavy analogs of A and they can be 
obtained from Eqs. (3) and (4) as previously. 

III. NUMERICAL SIMULATION 
A. Actions 

The gauge action as well as the heavy quark NRQCD 
action used for this work are described in detail in 
ReL The gauge action is tadpole improved and the 
leading classical error is quartic in lattice spacing. The 
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FIG. 1: Effective Mass M{t) versus t for singly heavy E-like baryons for different combinations of light and heavy quark mass 
(denoted by hopping parameter k and bare mass m, respectively) . Open symbols are for calculations with a correlation function 
with local source and sink, filled symbols are for local source and smeared sink. 



Hamiltonian corresponding to the NRQCD action is com- 
plete to 0{1/M^) in the classical continuum limit. For 
light quarks we use a Dirac- Wilson action of the D234 
type Q which has been used previously and detailed 
in Refs. j^, ||. Its leading classical errors are cubic in 
lattice spacing. All these actions are summarized in the 
appendix. 



B. Simulation Details 

This work is done with two sets of quenched gauge 
configurations (at /? = 2.1 and 2.3) on anisotropic lattices 
with a bare aspect ratio Us/at = 2, where the spatial 
lattice spacing varies from about 0.22 to 0.15 fm. 

The renormalized anisotropy is obtained from 



a, 



s. = — 

at atV{r2 



V{r2) - a,y(ri) 



(7) 



where V{r) is the potential between a static quark- 
antiquark pair with separation r, and is extracted from 
an exponential fit to a sequence of Wilson loops. In 
the numerator of Eq. (^ the sequence of Wilson loops 
extends in a coarsely-spaced direction, and in the de- 
nominator the sequence extends in the finely-spaced di- 
rection. The separation r may be along a lattice axis 
or off-axis, and various possibilities were included in the 



calculation. However, the separation r never includes the 
finely-spaced direction, so that Vir) itself is independent 
of as I at- It is convenient to avoid using the largest val- 
ues of r, where the exponential fit becomes noisier and 
the periodicity of the lattice can affect a determination 
of the anisotropy. Our results are 



f - ^ - i 1-96(2), for 13 = 2.1 
at~\ 1.99(3), for /3 = 2.3. 



(8) 



We used fixed time boundaries to construct quark 
propagators, and gauge fields were generated using a 
pseudo-heat-bath Monte Carlo algorithm with 400 (/3 = 
2.1) to 800 {(3 = 2.3) sweeps between saved configura- 
tions. For (3 — 2.1, we use 720 configurations and for 
/3 = 2.3 the number of configurations is 442. Two sets 
of bare masses are used for each heavy quark while four 
sets of hopping parameters are used for the light one. 
Bare masses for heavy quarks are chosen to surround the 
physical value so that an interpolation can be used. For 
example, at /3 — 2.1, the charm mass is in between 1.2 
and 1.5 and the bottom mass is in between 5.0 and 6.0. 
The charm mass is fixed by setting the rjc mass to its ex- 
perimental value, whereas the B° mass is used to fix the 
bottom mass. The hopping parameter corresponding to 
the strange quark is fixed from the Dg meson mass. The 
temporal lattice spacing and correspondingly the scale 
is fixed by setting the p-meson mass to its experimental 
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FIG. 2: Effective Mass M{t) versus t for doubly heavy E-like baryons for different combinations of light and heavy quark 
mass (denoted by hopping parameter k and bare mass m, respectively). Open symbols are for calculations with a correlation 
function with local source and sink, filled symbols are for local source and smeared sink. 



TABLE III: Summary of lattice parameters. The quantity 
a^^ is the inverse of the temporal lattice spacing while Us 
and ut are the tadpole improvement factors for spatial and 
temporal links respectively. 
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TABLE IV: Hopping parameters and bare masses. Four k 
values were used in simulations at each (3. Ks is the hopping 
parameter for the strange quark, and c and b are the charmed 
and bottom bare masses respectively. 
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value. Summaries of lattice parameters as well as hop- 
ping parameters for heavy and light fields are given in 
Tables III and IV respectively. 

Correlation functions are calculated using interpolat- 
ing operators in local form at both source and sink. In ad- 
dition to that we use a gauge invariant smearing for quark 
propagators at the sink using the smearing function from 
Eq. (13) of Ref. [Q. These local and sink-smeared corre- 
lators are fitted simultaneously to obtain hadron masses. 
The required correlations among different quantities are 
taken into account by covariant matrices obtained from 
singular value decomposition, and the statistical error is 
estimated from bootstrapping the fitting procedure. As 
in Ref. local correlators are fitted with two exponen- 
tial functions (Aexp (—mix) -l-Bexp (—77122;)), while the 
sink-smeared correlation function is fitted with a single 
exponential (Cexp (— mix)). The mass parameter for 



the sink-smeared fit is constrained to be the same as the 
lowest mass of the fit to the local correlator. The time 
window for the fit is chosen in a way such that the end- 
ing time is large and the fit is stable under variation of 
both starting and ending time by a few time steps. Light 
quark extrapolation is done by extrapolating the hadron 
masses extracted at four light quark masses with the form 
Co -l- C2m^ -|- cam^, where 171^^ is the pion mass. In most 
of the cases the cubic (m^) contributions are small and 
they are included only to get systematic errors. 

Figs. 1-3 show some representative examples of our 
simulation results. We plot the effective mass for dif- 
ferent heavy baryons versus time t, where the effective 
mass is defined to be M{t) = ln{g{t)/g{t + 1)) with g{t) 
being the zero-momentum time correlation function of 
baryon fields. Open symbols in these figures are for cal- 
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FIG. 3; Effective Mass M{t) versus t for doubly heavy E*-like baryons for different combinations of light and heavy quark 
mass (denoted by hopping parameter k and bare mass m, respectively). Open symbols are for calculations with a correlation 
function with local source and sink, filled symbols are for local source and smeared sink. 



culations with a correlation function with local source 
and sink, filled symbols are for local source and smeared 
sink. There is good agreement between local and smeared 
results at large times. 

It should be noted that the actual fits to determine the 
masses are performed directly with the correlation func- 
tions, and not on the effective masses plotted in Figs. 1-3, 
but the plots provide an indication of the quality of our 
data. Although the sink-smeared results appear to be 
somewhat noisy they are quite helpful in constraining 
the two-exponential fit of the local correlation function. 



C. Mass Extraction 

The kinetic mass of a nonrelativistic state can be ex- 
tracted from the usual NRQCD relation pT[ 

27r2 

2 

which is derived from E — tttt — . Here N ~ La^, with L 
being the lattice size and the spatial lattice spacing. 
^ = atlas is the anisotropy whereas Eq and Ep are simu- 
lation energies corresponding to the ground state and the 
state with momentum p — respectively. Mass dif- 
ferences between two states [Hi and H2) with the same 



heavy quark can be obtained by taking the difference of 
their zero momentum simulation energies: 

Mh.^Mh, = El,^{0)-El^{0), (10) 

which follows from the lattice NRQCD expression for the 
hadron mass 

AI„ - Es^„^{0)+ZMQ-Esh^ft, (11) 

where Egim is the simulation energy at zero momentum 
and the last two terms represent the renormalized heavy 
quark mass. The bare quark mass AIq has both a mul- 
tiplicative (Z) and additive (Eshift) renormalization[p^ 
which should be independent of hadronic state. A more 
precise result is obtained for heavy hadron masses with 
the heaviest light quark (k = 0.229) rather than with a 
lighter fight quark {k = 0.233 and higher). 

Moreover, mass differences (Eq. 10) can be calculated 
more precisely than masses (Eq. 9). Therefore, for exam- 
ple, one can calculate a meson mass from the relation 

M{qi, Q) = M{qh,H) - AM =. M{qH,H) - AE, (12) 

where 

AM^AE = Eiqh,Q)-E{qi,Q). (13) 

Here qh and qi denote the heaviest light quark and a 
lighter one respectively, and M{qh,H) is extracted by 
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TABLE V: Results for meson and baryon masses and mass 
splittings (in MeV) calculated using the NRQCD action for 
charm (c) and bottom (b) quarks. The first error is statis- 
tical while the second error comprises systematic errors due 
to scale, time window and anisotropy. Rows are separated 
into mesons, singly heavy charmed baryons, doubly heavy 
charmed baryons, singly heavy bottom baryons and doubly 
heavy bottom baryons, respectively. 
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using Eq. (9). Eq. (13) is valid as long as Z in Eq. (11) 
is the same i.e., both states consist of the same heavy 
quark Q. 

Similarly, masses of singly and doubly heavy baryons 
can be extracted from meson masses by using 
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where 



Esh - E{qh,Q)~ E{qiq2,Q), 
Edh = E{QQ)-E{qi,QQ). 
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TABLE VI: Results for charmed baryon masses and mass dif- 
ferences (in MeV) compared to experimental values. The first 
row of lattice results (taken from |^) were calculated using 
a relativistic action of the D234 type for the charmed quark 
while for the second row (this work), the NRQCD action was 
used. 
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For example, the mass is extracted by taking its dif- 
ference (at each k) with the D{B^) mass (m) at k = 0.229 
and then subtracting that from to. Masses extracted by 
using Eq. (9) and Eqs. (12-17) are consistent with each 
other. However, errors in the second method are smaller 
than the previous one. 



IV. RESULTS 

The mass spectrum and spin splittings of heavy quark 
baryons have been computed on an anisotropic lattice 
using the NRQCD heavy quark action. Results are sum- 
marized in Table V, where the first error is the statistical 
error obtained from a bootstrap analysis with a boot- 
strap sample size equal to the configuration sample size. 
The second error is an overall systematic error due to 
scale and anisotropy uncertainties, the uncertainty due 
to choosing a time window, the light quark extrapolation 
error and the strange quark mass uncertainty. Mesons, 
singly heavy baryons and doubly heavy baryons are sep- 
arated into different groups by horizontal lines. In Table 
VI we have compared our results with those obtained by 
using a relativistic (D234) heavy quark action and ex- 
perimental numbers (where available). One can notice 
that the NRQCD results and D234 results are consistent 
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FIG. 4; Spin splittings in the meson sector, plotted as 
{Mv — Mp)M versus M where M is the average vector and 
pseudoscalar meson mass [Mv + Mp)/2. "Lattice literature" 
results are from Refs. |0, M, 0| while D234 results are from 
Ref §. 



FIG. 5; Spin splittings in the baryon sector plotted as 
[M-j,/2 — Mi/2)M versus M where M is the average vector 
and pseudoscalar meson mass {Mv + Mp)/2. "Lattice liter- 
ature" results are from Ref. |lq] while D234 results are from 
Ref |8|. The solid line is a fit to the experimental data 



with each other. Results are also consistent with a pre- 
vious NRQCD calculation Q. As in Ref. it is found 
that the suppression of spin splittings is not present in the 
baryon sector, although such a suppression is known to 
be characteristic of the heavy meson sector. One can also 
notice that the spin splittings for doubly heavy baryons 
are as large as their singly heavy counterparts. 



V. DISCUSSION AND SUMMARY 

In order to put the results of the present calculation 
into perspective it is useful to consider spin splittings 
over the whole range of available quark masses. We start 
with mesons where it has been known for a long time 
that the squared mass difference My — Mp for vector 
and pseudoscalar mesons is approximately constant for 
all mesons of the form Qq, where q is up or down and Q 
is any light or heavy flavor. This relation is illustrated 
in Fig. 4 for the mass pairs (p, tt), {K*,K), {D* , D) and 
{B* , B). Also shown are the results of lattice simulations 
including the present work. The tendency for quenched 
lattice QCD to underestimate the spin splittings relative 
to experimental values is clearly visible. 

In Ref. we showed that it is useful to consider the 
behavior of the spin splittings in the baryon sector as 
a function of quark mass also in terms of the mesonic 
average mass {My + Mp)/2. The results for the baryon 
pairs (A, TV), (E*,E), (S*,Sc) and (SJ,Sb) are shown in 
Fig. 5. It is a remarkable empirical fact that the baryon 
spin splitting scales almost exactly like the inverse of the 
average meson mass. The implication is that the ratio 
of meson to baryon spin splittings is almost constant. 
This was discussed in Ref. [B and was to some extent 



anticipated by Lipkin from the point of view of the 
quark model (see also Lipkin and O'Donnell |]l9| ). 

The results of quenched lattice calculations are also 
shown in Fig. 5. The suppression of spin splittings rel- 
ative to experiment, visible for mesons, is not seen for 
baryons. The results of the present lattice NRQCD cal- 
culation support this conclusion in the charm and bottom 
sectors. It is clear that a definitive measurement of Ef, 
and masses would be highly desirable to extend the 
experimental comparison to larger mass values. 

From the point of view of lattice NRQCD our re- 
sults present an interesting challenge. As is well known, 
the spin splittings of both charmonium |Q and heavy- 
light mesons |l^ are clearly underestimated by lat- 
tice NRQCD simulations. Up to now, these simulations 
have used couplings modified only by mean-field tadpole 
factors. It is tempting to speculate that there are addi- 
tional large perturbative corrections to these couplings. 
In particular it might seem that a correction to the quark 
coupling with the chromomagnetic field (the C4-term in 
Eq. (A5)) has the potential to cure the spin splitting de- 
ficiency for both charmonium and heavy-light mesons. 
However one has to be cautious in wishing for such a 
cure as it would upset the already reasonable values for 
spin splittings in the baryon sector. 

To summarize, we have calculated the masses of 
baryons containing one or two heavy quarks using 
quenched lattice QCD. NRQCD is used to describe charm 
and bottom quarks. In the charm sector the results of 
this work are compatible with those obtained previously 
where a Dirac- Wilson action of the D234 type was used 
for the heavy quark. No suppression of the spin split- 
tings observed in lattice NRQCD simulations of heavy- 
light mesons is seen in the heavy baryon sector. 



8 



This and our previous work |^ leave a number of diffi- 
cult open questions. One would like to be able to improve 
the lattice calculations of baryons to reduce the uncer- 
tainties to the same level achievable in mesons. Also 
how (and whether) the addition of dynamical quarks to 
the simulations will solve the dilemma of spin splittings 
has yet to be understood. A phenomenological issue is 
to understand the remarkable constancy in the meson to 
baryon spin splitting ratio over the whole available quark 
mass range. On the experimental side it will be a signif- 
icant challenge to extend baryon mass measurements in 
the bottom and doubly heavy sectors. 
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APPENDIX A: DETAILS OF ACTIONS 

1. NRQCD Action 

The heavy quark action is nonrelativistic and is dis- 
cretized to give the following Green's function propaga- 
tion: 



/ atHB\ A atHA \ ' 



ut 



2n 



1-^|G.,(A1) 



The nonrelativistic Hamiltonian is complete to 0{1/M^) 
in the classical continuum limit : 



H = Ho + SH, 
-A(2) 



Hn = 



2M 



(A2) 
(A3) 



C4 ^ a?A(4) 



u^2M ^ 24M 



(A5) 



u^ul SAP ^ 

C3 9 



7,(T • (A X E - E X A) 

uWf 8M2 ^ ' 



C6 , (A6) 



-ci 



16<M2 ' 



8A/3 ui 8Af3 



8M3 





E X E 








B2\ 







ali^' ,,,, 

"^"192^2^2^/3 • ^^^^ 



Here a tilde signifies discretization errors have been re- 
moved. In particular. 



(A8) 
(A9) 



6^2 



E^ = 




= 


2^ijkFjk 




If,a^) 







The various spatial lattice derivatives are defined as 
follows: 

a,A,G(x) = -^[Ui{x)G{x + i) 

-Ul{x-i)G{x-i)], (All) 
asA[+^G{x) = ^^^G{x + i)-G{x), (A12) 

a,A';-^G{x) = G(x) - ' G{x - i), (A13) 

Us 

alA^^G{x) = ^G(x + £) - 2G(x) 



U^x - i) 



G{x-i), (A14) 



A, = A, - ^a|+)a,A^\ (A15) 



A(2) = ^Af , (A16) 



A(2) ^ A(2)_^A(4), (A17) 
12 ' ^ ' 

A(4) = ^(Af))'. (A18) 



2. Gauge Field Action 

The leading classical errors of the gauge field action 
are quartic in lattice spacing. The action is 



SaiU) 



5/3 

T 



i^2:(i-iR,w„) 



pt 
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20uii 



- -RcTrC/rst 



where anisotropic ratio ^ = a^/at and /3 is the lattice 

gauge field coupling constant. 

ps : spatial plaquettes 

rs : spatial planar 1x2 rectangles, 

pt : plaquettes in the temporal-spatial plane, 

rst(rts) : rectangles with the long side in a spa- 

tial(temporal) direction. 



3. Light Quark Action 

For light quarks, we used a D234 action |^ |l^ with 
parameters set to their tadpole-improved classical values. 
Its leading classical errors are cubic in lattice spacing and 
the action can be written as 



4k ^ 



1 



4k 



E 



1 



Du{x) 



1 



1 



2^ 



+ 7r4^ Y. i'ix)ao^Fo^ix)iP{x) - V ^(x)^(x), (A20) 



(A19) where 



Du{x) = '4>{x){l - ^-f^)U,{x)ilj{x + i) 

+i;ix + i)il+^-f,)uUxMx), (A21) 
Duix) = ^ljix){l-j4)Ui{x)ij{x + i) 

+^(x + i)(l + 74)C/i(a;)V'(a;), (A22) 
D2^{x) = il>{x){l'^j,)U,{x)U,{x + i)^j{x + 2l) 

+ 4'{x + 2i){l + ^-fi)U^^{x + i){7/(x)^(x), (A23) 
D2t{x) = 4>{x){l - j4)U4{x)U4{x + i)i:{x + 2i) 

+ {l^{x + 2i){l + j4)Ul{x + i)ul{x)il^{x), (A24) 

gF^Ax) - ^^{n^.{x)-nl,{x)) 



-Im(Trl]^,(x)), 



(A25) 



- —[U^{x)UA^ + fi)Ul{x + i))Ul{x) 



+ Uu{x)Ul{x - fL + v)Ul{x - il)U^,{x - fi) 
+ Uj^{x — (i)Ul{x — (i — v)U^i{x — A — v)Uu{x — v) 
+Ul{x - v)U^{x - v)U,{x + fi- 0)Ulix)] . (A26) 



x.i<j 
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